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CN ' Abstract. A new method is presented for modeling the transformation 

! between two polarimetric pulse profiles in the Fourier domain. In practice, 

<D • one is a well-determined standard with high signal-to-noise ratio and the 

other is an observation that is to be fitted to the standard. From this fit, 
' both the longitudinal shift and the polarimetric transformation between 

CN , the two profiles are determined. Arrival time estimates derived from the 

best-fit longitudinal shift are shown to exhibit greater precision than those 
, derived from the total intensity profile alone. In addition, the polarimetric 

' transformation obtained through this method may be used to completely 

. calibrate the instrumental response in observations of other sources. 
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O ■ !• Polarimetric Pulse Profile Modeling 

Q^. In radio astronomy, a pulsar's mean polarimetric pulse profile is measured by 

averaging the observed Stokes parameters as a function of pulse longitude. By 
integrating many well-calibrated pulse profiles, a standard profile with high 
I signal-to-noise ratio (SNR) may be formed and used as a template to which 

■ individual observations are fit. For example, in Appendix A of Taylor (1992), 
^ I a method is presented for modeling the relationship between standard and ob- 

■ served total intensity profiles in the Fourier domain. In the current treatment, 
rS I the scalar equation that relates two total intensity profiles is replaced by an 

■ analogous matrix equation, which is expressed using the Jones calculus. The 
polarization of the electromagnetic field is described by the coherency matrix, 
p = (^o (To + S • <t)/2, where 5o is the total intensity, S = {Si, S2, S3) is the 
Stokes polarization vector, ctq is the 2x2 identity matrix, and cr = (cri, <T2, cr^) 
are the Pauli spin matrices (Britton 2000). Under a linear transformation of the 
electric field vector as represented by the Jones matrix, J, the coherency matrix 
is subjected to the congruence transformation, p' = JpJ^ (Hamaker 2000). 

Let the coherency matrices, p'{(j)n), represent the observed polarization as a 
function of discrete pulse longitude, (pm where < n < and N is the number 
of pulse longitude intervals. Each observed polarimetric profile is related to the 
standard, Poifpn), by the matrix expression, 

p'iM = PdC + PN(0n) + Jpo('^n " '/')J^ (1) 

where p^Q is the DC offset between the two profiles, represents the system 
noise, J is the polarimetric transformation and ip is the longitudinal shift between 
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the two profiles. The Jones matrix, J, is analogous to the gain factor, b, in 
equation (1) of Taylor (1992). However, as J has seven non-degenerate degrees 
of freedom, the matrix formulation introduces six additional free parameters. 
The discrete Fourier transform (DFT) of equation (1) yields 

p'(t'm) = PN(f^m) + Jpo{i^m)J^ exp{-i27rmcp) , (2) 

where is the discrete pulse frequency. Given the measured Stokes parameters, 
S'k{4'n), and their DFTs, S'^(fm), the best-fit model parameters will minimize 
the objective merit function, 



= 1 k=0 



where <jjt is calculated from the noise power and Tr is the matrix trace. As in 
van Straten (2004), the partial derivatives of equation (3) are computed with 
respect to both ip and the seven parameters that determine J. The Levenberg- 
Marquardt method is then applied to find the parameters that minimize x^- 



2. Applications 

Dual-polarization observations of PSR J0437— 4715 were made using the Parkes 
Multibeam receiver and CPSR-II, the 128 MHz baseband recording and real-time 
processing system at the Parkes Observatory. The data were observed during 
seven separate sessions between 5 June and 21 September 2003, calibrated using 
the method described in van Straten (2004), and integrated to produce the 
polarimetric standard shown in Figure 1. 

2.1. High-Precision Timing 

In any high-precision pulsar timing experiment, the confidence limits placed on 
the derived physical parameters of interest are proportional to the precision with 
which pulse time-of-arrival (TOA) estimates can be made. Aside from typical 
observational constraints such as system temperature, instrumental bandwidth, 
and allocated time, TOA precision also fundamentally depends upon the physical 
properties of the pulsar, including its flux density, pulse period, and the shape 
of its mean pulse profile. When fully resolved, sharp features in the mean pulse 
profile generate additional power in the high frequency components of its Fourier 
transform. As higher frequencies contribute stronger constraints on the linear 
phase gradient in the last term of equation (2), sharp profile features translate 
into greater arrival time precision. 

This important property may be exploited in order to significantly improve 
the precision of arrival time estimates derived from full polarimetric data. As 
noted by Kramer et al. (1999), the mean profiles of Stokes Q, U, and V may 
contain much sharper features than that of Stokes /, especially when the pulsar 
exhibits transitions between orthogonally polarized modes. Although the po- 
larized flux density is lower than the total intensity, these transitions lead to a 
greater SNR in the high frequency components of the polarized power spectra, 
as shown in Figure 2. 
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Figure 1. Mean polarization of PSR J0437— 4715, plotted using po- 
lar coordinates: orientation, 9; ellipticity, e; and polarized intensity, 
5" = (plotted in red below the total intensity). Flux densities are 
normalized by ao, the r.m.s. of the off-pulse total intensity phase bins. 
Data were integrated over a 64 MHz band centered at 1341 MHz for 
approximately 40.6 hours. 



To demonstrate the potential for improved arrival time precision, TOA esti- 
mates spanning over 100 days were derived from the 490 five-minute integrations 
used to produce the standard plotted in Figure 1. The arrival times derived from 
the mean total intensity profile have a post-fit residual r.m.s. of 198 ns. By mod- 
eling the polarimetric pulse profile as described in Section 1, the resulting arrival 
times have a post-fit residual r.m.s. of only 146 ns, an improvement in precision 
of approximately 36%. 



2.2. Instrumental Calibration 

Based on the assumption that the mean polarimetric pulse profile does not vary 
significantly with time, the standard profile and modeling method may also be 
used to determine the polarimetric response of the observatory instrumentation 
at other epochs. As a demonstration, a single, uncalibrated, five-minute inte- 
gration of PSR J0437— 4715 was fitted to the polarimetric standard shown in 
Figure 1. The model was solved independently in each of the 128 frequency 
channels, producing the instrumental parameters shown with their formal stan- 
dard deviations in Figure 3. In each 500 kHz channel, it is possible to estimate 
the ellipticities and orientations of the feed receptors with an uncertainty of 
only one degree. This unique model of the instrumental response may be used 
to calibrate observations of other point sources. 
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Figure 2. Power spectra of the standard profile plotted in Figure 1. 
Stokes /, Q, U , and V are plotted in black, red, green, and blue. 
Above 40 times the pulse frequency, both Stokes Q and V exhibit a 
SNR greater than or equal to that of the total intensity. 

3. Conclusion 

When compared with the scalar equation used to model the relationship between 
total intensity profiles, the matrix equation presented in Section 1 quadruples the 
number of observational constraints while introducing only six additional free 
parameters. By completely utilizing all of the information available in mean 
polarimetric pulse profiles, arrival time estimates may be obtained with greater 
precision than those derived from the total intensity profile alone. In addition, 
the modeling method may be used to uniquely determine the polarimetric re- 
sponse of the observatory instrumentation using only a short observation of a 
well-known source. 
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Figure 3. Best-fit instrumental parameters as a function of ob- 
serving frequency. A single, uncalibrated, five-minute integration of 
PSR J0437— 4715 was fitted to the polarimetric standard. From top 
to bottom are plotted the orientations, 9^, and ellipticities, e^, of the 
feed receptors, the differential phase, (j), the differential gain, 7, and 
the absolute gain, G, specified in units of an intermediate reference 
voltage that can later be calibrated to produce absolute flux estimates. 
In the top two panels, black and red correspond to receptors and 1, 
respectively. 



